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The  Aggregate  Behavior  of  Branch  Points 
-  Branch  Point  Density  as  a  Characteristic 
of  an  Atmospheric  Turbulence  Simulator 

Denis  \V.  Oesch“,  Darryl  J.  Sanchez",  Carolyn  M.  Tewksbury-Christle",  Patrick  R.  Kelly" 
"Starfire  Optical  Range,  Air  Force  Research  Labs,  Kirtland  Air  Force  Base,  New  Mexico,  USA 


ABSTRACT 

The  Atmospheric  Turbulence  Simulator  used  in  testing  in  the  Atmospheric  Simulation  and  Adaptive-optic  Lab¬ 
oratory  Test-bed  at  Starfire  Optical  Range  is  configured  based  on  three  characteristics;  Fiied’s  parameter,  tq, 
the  Rytov  number,  and  the  Greenwood  Frequency,  fc-  All  three  may  be  estimated  from  open  loop  data 
as  a  means  of  verifying  the  simulated  turbulence  conditions  for  a  given  test  configuration.  However,  unlike  tq 
and  fc,  the  Rytov  number  isn’t  directly  calculated.  Instead  the  scintillation  index  is  estimated  from  intensity 
measurements.  At  low  Rytov  values,  (<  0.3  —  0.4),  this  measurement  can  approximate  the  Rytov  number,  how¬ 
ever  beyond  a  Rytov  of  0.4  this  parameter  becomes  saturated.  Branch  Points  begin  to  appear  after  the  Rytov 
value  exceeds  0.1.  In  this  work  the  behavior  of  the  branch  point  density  is  examined  to  determine  its  viability 
as  another  parameter  for  calibration  our  turbulence  simulator. 

Keywords:  branch  points,  density,  atmosphere,  turbulence,  adaptive  optics 

1.  INTRODUCTION 

To  calibrate  our  experimental  conditions  in  the  Atmospheric  Simulation  and  Adaptive-optic  Laboratory  Testbed 
(ASALT)  at  Starfire  Optical  Range,  it  is  standard  practice  to  collect  open  loop  data  during  experiments  for  each 
Atmospheric  Turbulence  Simulator  (ATS)'  configuration  under  which  closed  loop  data  is  collected.  This  open 
loop  data  provides  an  independent  verification  that  the  tests  were  conducted  under  the  prescribed  conditions. 
Processing  of  this  data  provides  measurements  of  Fried’s  parameter,  ro,  the  Greenwood  Frequency,  fc,  and  the 
scintillation  index. 

The  scintillation  index  is  approximately  equivalent  to  the  Rytov  parameter  for  values  of  the  Rytov  <  0.3  — 0.4. 
It  is  well  known  from  experimental  data  that  at  Rytov  0.1  branch  points  begin  to  appear  in  the  data.  It  was 
suggested  that  branch  point  density  might  be  used  as  an  alternative  to  the  scintillation  index  for  calibrating  the 
configuration  of  the  ATS.  The  purpo.se  of  this  work  was  to  examine  the  functional  dependencies  of  the  branch 
point  density  on  other  system  parameters  to  determine  it’s  u.se  as  a  calibration  parameter. 

2.  BACKGROUND 

The  ASALT  laboratory  is  uniquely  suited  to  studying  the  saturation  regime,  Rytov  >  0.4.  The  laboratories 
have  3  full  bench  top  adaptive  optic  systems  each  with  its  own  ATS  and  multiple  deformable  mirror  (DM)  - 
wavefront  sensor  (WFS)  pairs,  all  operating  at  1.55  /im.  For  this  work  the  WFS  was  a  temporal  SRI  with  256x256 
pixel  resolution  over  a  simulated  1.5  m  aperture.  The  ATS  consists  of  two  6  inch  phase  screens  etched  with  a 
Kolmogorov  structure  function  set  in  converging  beam  space  providing  for  a  range  of  turbulence  conditions.  By 
varying  the  position  of  each  phase  wheel  along  a  rail  provides  a  method  for  adjusting  tq  between  roughly  4  and 
30  cm  while  also  varying  Rytov  from  0  up  to  roughly  2.4.  Stepper  motors  control  the  rate  of  rotation  of  the 
phase  wheels  and  allow  for  selectable  Greenwood  Frequencies,  fc-  The  wheels  are  fitted  with  a  hardware  home 
that  provides  repeatability  of  turbulence  conditions  over  multiple  tests.  To  this  standard  ATS  set-up  an  optical 
trombone  was  added  between  the  low  altitude  phase  wheel  and  the  system  pupil  to  provide  for  a  variable  free 
space  propagation. 

Further  author  information:  (Send  correspondence  to  Denis  W  Oesch) 

Ei-mail:  denis.w.oesch@saic.com 
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The  optical  trombone  allows  partial  decoupling  of  the  ro  and  Rytov  parameters.  Thus  providing  a  means 
of  varying  the  distance  to  the  turbulence  layer  and  hence  the  Rytov  parameter,  while  maintaining  a  constant 
turbulence  strength.  This  made  isolating  the  contributions  to  the  distribution  of  branch  points  seen  in  the 
pupil  plane  from  turbulence  strength  and  distance  independently.  Leading  to  an  understanding  of  how  the 
configuration  of  the  phase  wheels  relates  to  the  measured  branch  point  densities. 

There  has  been  significant  work  on  branch  points  in  atmospheric  turbulence  and  their  impact  on  adaptive 
optics  systems^"*^  and  even  one  specifically  on  branch  point  density.'  *  However  these  hasn’t  been  an  experimental 
test  of  this  t>T>e  that  we  have  found  in  the  literature  which  allows  for  the  isolation  of  the  turbulence  strength 
and  the  propagation  distance  as  has  been  done  in  this  work. 


3.  PROBLEM  FORMULATION 

The  Rytov  approximation  forms  the  theoretical  basis  of  adaptive  optics.*"  With  this  approximation  and  a 
Kolmogorov  structure  function,  it  has  been  shown*"’  *®  that  four  moments  of  the  atmospheric  structure  function 
determine  AO  performance,  the  zeroth  altitude  moment,  ro,  the  (5/6)*^*  altitude  moment,  (tJ,  the  (5/3)*’^ 

altitude  moment,  Oq,  and  the  (5/3)*^^  velocity  moment,  fc-  We  generate  these  parameters  using  our  atmospheric 
turbulence  simulator  (ATS). 

Our  ATS  consists  of  two  turbulence  layers  each  with  adjustable  strengths  and  altitudes.  In  its  standard 
configuration,  ro  and  <t^  are  coupled  because  both  are  functions  of  ATS  phase  wheel  position  (and  phase  wheel 
position  Fresnel  scales  to  altitude).  To  partially  mitigate  the  coupling  of  these  parameters,  an  optical  trombone 
was  added  in  collimated  space  after  the  second  (low  altitude)  phase  wheel.  This  allows  for  an  increase  in 
propagation  distance  without  changing  the  atmospheric  structure  function  and  this  allows  variation  of  without 
changing  ro,  i.c.  the  altitude  of  the  phase  screens  can  be  displaced  by  AL  causing  — ♦  2)^^  +  AL  and 
^high  “high  ■*"  "'*'**  low  ®***^  high  denoting  the  low  and  high  altitude  structure  functions, 

respectively.  Then, 

^n(^)  =  “  ^low)  ~  ^high) 

which  causes  the  no  trombone  version  of  the  Rytov  parameter, 

al=0.291kl^'^  rc2(2)(v)"/«d2 

Jo 


to  go  to 

(^x}with  trombone 


=  C2(2)(72)"''"d2 


=  (<?;j,„,(%w)(-«(2low  +  At))"'"  +  ,,,gh(*higl,)(T'(%Bl,  +  At))*'") 

^  ^‘^xJno  trombone 


but  the  no  trombone  version  of  the  Fried  parameter. 


(2) 


To 


=  ^A23k^J^ 


CUz)dz 


-3/r, 


goes  to 


{'•o}with  trombone  =  (0-423A-2  C?,(z)  dz) 


-3/5 


-  (0.423fcg  (C2  |o^v(2l 
~  {*’o}no  trombone 


)  +  C 


n.high 


(^high))) 


-3/5 


(3) 


where  7  =  1  for  a  plane  wave,  ko  the  wave  number,  L  either  the  distance  to  the  guide  star  or  the  top  of 
the  atmosphere,  and  the  atmospheric  structure  function.  This  is  equivalent  to  the  electromagnetic  field 
encountering  turbulence  followed  by  a  free  space  propagation  of  variable  length.  The  separation  of  the  turbulence 
layer  and  the  pupil  can  be  adjusted  which  changes  the  Rytov  parameter  and,  in-turn,  alters  the  number  of  branch 
points  ill  the  pupil. 
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3.1  Definition  of  a  Branch  Point 

A  branch  point  is  a  zero  of  an  analytic  function.  Following  FVied,^’®  let  //)  be  the  scalar  function  representing 
an  optical  field  and  let  u{x,  y)  and  u(j:,  y)  be  the  real  and  imaginary  parts  of  that  function,  respectively,  then 

w{x,  y)  =  u(t,  y)  -f  >a(x,  y).  (4) 

where  the  amplitude,  A(x,i/),  and  phase,  0{x,y),  of  this  field  will  be  given  by 

A{x,y)  -  \/u^{x,y)  +v^{x,y)  (5) 

and 

d(x.  y)  =  ar!i(w{x,  y))  +  2ffK.  (6) 

Where  k  =  k(x,  y)  are  two-dimensional  step  functions  of  height  25r. 

Prom  the  theory  of  analytic  functions,  a  branch  point  at  (xp,  2/p)  is  equivalent  to  both 

iwCxp,  yp)  =  0  (7) 

and 

f  dF  V^(x'.y')  =  ±27r  (8) 

Jc 

for  C  a  closed  curve  encircling  the  one  branch  point  at  (xp,?/p).  and  V4,  the  gradient  operator  acting  on  <{).  Also, 
for  discussional  clarity,  we  call  sign  ^  dl'  ■  V,i{x',y')  =  ±27r^  the  polarity  of  the  branch  point. 

3.2  Identification  of  a  Branch  Point’s  Position  and  Polarity  in  Experimental  Data 

Experimentally,  the  analytic  function  is  sampled,  such  that  the  measurement  is  an  average  over  a  detector  pixel. 
Therefore,  the  theoretical  expressions  given  in  F.quations  7  and  8  must  be  modified  such  that  given  a  branch 
point  at  (xp,2/p)  and  pixels  centered  at  Pa  :=  (x^.i/^),  Pq  :=  (x^  +  AX,yA),  Pc  ■=  i^A  +  ^X,yA  -  AT),  and 
Pd  :=  {xa-^Va  —  AF)  such  that  x,4  <  Xp  <  x^  -i-  AX  and  yA  <  Vp  <  Va  +  AF  then 

7^  0  (9) 

where  ic,,,  is  the  measured  value  of  the  scalar  function,  w(xp,  yp).  Hence,  an  experimental  measurement  of  the 

zero  crossings  of  w{x,y)  cannot  be  used  to  determine  the  location  of  branch  points.  On  the  other  hand,  the 
circulation  in  phase  given  by  the  integral  in  Equation  8.  dU  ■  V^{x',  y').  when  sampled,  goes  to 

(<?(Pb)  -  <!>(P^))  +  {o{Pd)  -  <>{Pb))  +  (o(Pc)  -  o(Pd))  +  (d(P4)  -  <p(Pc))  (10) 

with  the  contour  C  given  by 

C'  :  (P.0  -  (Pb)  -  (Pd)  -  (Pc)  -  (P.4),  (11) 

(see  Figure  1).  Equation  10  still  measures  the  ±27r  circulation  in  phase  indicative  of  branch  points,  and  hence, 
this  expression  is  itsed  to  find  branch  points  in  experimental  data.  To  implement  this,  we  sum  the  gradients  of 
the  phase  around  a  closed  loop  in  the  2ir  modulo  phase,  with  all  loops  chosen  to  be  elementary  circulations,  i.e. 
a  2x2  region  of  the  27r  modulo  phase  was  used  for  each  loop.  If  the  sum  of  the  phase  gradients  is  ±2;r,  an  odd 
number  of  branch  points  exists  within  the  loop.  We  assume  that  there  is  only  one  branch  point  within  the  loop. 
We  then  assign  the  location  of  the  branch  point  as  the  intersection  of  the  four  pixels  with  its  polarity  given  by 
the  sign  of  the  circulation.  This  is  a  standard  technique  for  isolating  branch  points  from  phase  data.^ 

4.  BRANCH  POINT  DENSITY 

The  density  of  the  branch  points  should  increase  with  propagation  distance  and  provide  additional  information 
in  characterizing  the  profile  of  the  atmospheric  turbulence.  To  determine  the  relative  influence  of  the  distance 
tmd  strength  of  the  turbulence,  open  loop  data  was  compiled  with  the  temporal  SRI.  The  data  collected  from 
this  camera  was  processed  for  branch  point  density  in  a  range  of  configurations. 
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4.1  Data  and  Methodology 

I'or  these  test  sequences,  a  single  phase  wheel  was  used  in  conjunction  with  the  optical  troinhone  descrihed  in 
Section  ??.  In  each  data  run.  the  strength  of  the  turbulence  is  maintained  while  the  distance  to  the  turbulence 
layer  is  increiised.  This  incr(’iuse  takes  the  form  of  eleven  evenly  spaced  steps  over  a  roughly  9  kin  range,  providing 
a  test  region  from  0-15  km.  This  independently  varies  the  Rytov  jiarameter  from  0  -  1.2  as  shown  in  Figure  2. 
Between  each  data  run.  the  turbulence  strength,  which  in  our  test-bed  is  equivalent  to  varying  ro,  is  modified. 
In  this  way.  eleven  \-aluc*s  of  atmospheric  strength  were  interrogated,  ranging  from  G.52  to  16.60  cm.  See  Table  1. 
For  each  tromlrone  position  and  turbulence  strength  the  phase  wiieel  wtrs  rotated  at  a  constant  speed  for  200 
frames.  Data  wius  collcctcfl  for  each  configuration  producing  200  —  2.56.c256  images  of  the  phase  for  each  of  121 
turlnilence  conditions. 


Figure  1.  Sample  '2-  moditlo  phase  with  enlarged  portion  shotting  an  ej:ainple  of  an  elementanj  circulation. 


Figure  2.  Theoretical  Rytov  valne-.s  for  the  .single  phase  wheel  test  configurations  as  displayed  in  Table  1.  Each  curve  is 
ivpir.sentativc  of  the.  predicted  Rytov  number  for  a  series  of  Immbone  po.sitions,  0  —  ."lO  cm  through  a  constant  turbulence, 
strength.  Here  iv.prc.sented  as  a  constant  m  as  indicated  in  the  legend.  All  values  aiv  computed  for  a  simulated  1..')  meter 
ajierture. 


4.2  Measured  Density 

The  tlensity  of  tlic  l)raiidi  points  is  easily  cakulalecl  by  simply  api)lyiiig  the  techniques  in  Section  3.2  to  find  the 
branch  points,  then  averaging  the  number  of  branch  points  over  all  frames  collected  and  dividing  by  the  area  of 
the  aperture.  This  is  done  for  each  of  the  121  configurations  given  in  Table  1.  For  each  of  the  data  sets,  Figure  3 
plots  branch  point  density  along  the  vertical  axis  versus  propagation  distance  on  the  horizontal  axis.  Each  curve 
represents  a  different  turbulence  strength  denoted  by  the  legtmd.  \ha  Ecination  3.  turbulence  strength  is  given 


Branch  ='oint  Density 


Figure  3.  Branch  point  density  from  sinylc  phase  wheel  data  for  seler.ted  lurlmlence  strengths.  For  each  curve,  the  turbulence 
strength,  given  by  ro,  is  held  constant  while  the  propagation  distance  is  varied  u.sing  the  optical  trombone.  Raw  branch 
point  density  vs  propagation  distance. 


as  ro  in  centimeters  in  output  space  iussmning  a  l.Srn  pupil. 

As  can  be  seen,  the  branch  point  density  increases  with  propagation  distance.  It  also  increases  with  ttirbulence 
strength.  Of  the  two,  there  is  a  stronger  dependence  on  turbulence  strength  than  propagation  distance. 

4.3  Transforming  by  zo 

The  curves  presented  in  Figure  3  appear  to  shift  to  the  left  as  tire  .strength  of  the  turbulence  increases.  Implying 
that  there  may  be  some  minimum  distance,  cq  under  which  branch  points  simply  won't  form  for  a  given  turbulence 
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(103.0,  55.5) 

(3.0) 

4.2-12.6 

(7.27,0) 

10.3 

0.19  -  0.48 

7 

200 

7 

{0:5:50} 

(102.5,  55.5) 

(3.0) 

4.6-12.9 

(6.87,  0) 

9,7 

0.23  -  0.54 

8 

200 

7 

{0:5:50} 

(101.5,  55.5) 

(3.0) 

5.2-13.5 

(6.20,  0) 

8.8 

0.30  -  0.67 

9 

200 

7 

{0:5:50} 

(101.0,  55.5) 

(3.0) 

5.4-13.8 

(5.91,0) 

8.3 

0.34  -  0.73 

10 

200 

7 

{0:5:50} 

(99.5.55.5) 

(3.0) 

6.1-14.4 

(5.18,  0) 

7.3 

0.46  -  0.95 

11 

200 

7 

{0:5:50} 

(98.0.55.5) 

(3,0) 

G.G-15.0 

(•l.Gl.O) 

6.5 

O.GO  -  1.19 

Table  1.  Turbulence  parameters  used  for  measurement  of  branch  point  density  and  separation.  For  each  Ensemble  set,  200 


frames  of  data  was  collected  for  each  of  II  positions  of  the  optical  trombone,  ranging  m  5  cm  increments  from  0  to  50  cm. 


This  corresponds  to  ranges  of  propagation  distances  and  Rytov  parameters  for  each  Ensemble. 


0 


strength.  From  ex|jerimeiilal  results,  the  first  l)rain:h  points  form  wlien  a;  0.1.  Ev'ahiating  Kquation  2.  for  a 
single  turbulence  layer  produces  the  relationship 


0.563 1( 


b.88^,  -5/6  -.V.l  0/6 

2dH  ■ 


(12) 


Setting  ftj  =  0.1  and  solving  for  /..  allows  Equation  12  to  he  transformed  into  an  e(|uivalent  distance,  ;o 

;o  =  O.O-MSA-or,-;.  (13) 


which  can  he  considered  tin?  theoretical  value  at  whicli  branch  points  begin  to  a|)pear.  Then  Equation  13  can 
be  used  to  calculate  zq  for  the  11  turbuhmce  cases,  'I'hen  the  proiragation  distance  is  shifKxi  by  cq-  To  see  the 
magnitude  of  this  effect,  the  data  in  Figure  3  is  repKitted  in  Figure  I  with  the  idol’s  domain  recast  into  c  -  Zq- 


Branch  Point  Density 


Figtire  4.  Dmnch  point  density  from  single  phase  wheel  data  for  selected  turbulence  strengths.  For  each  curi'e.  the  turbulence 
strength,  given  by  ro,  is  held  constant  while  the  propagation  distance  is  varied  using  the  optical  trombone.  Branch  point 
density  vs  propagation  distance  with  plotted  versus  sealed  propagation  distance  (z  —  zqJ. 


This  fortn  clearly  demonstrates  that  branch  point  density  hits  a  stronger  de])endcnce  on  turbulence  strength 
than  propagation  distance. 

4.4  Normalizing  by  ro 

The  i.so-ro  density  curves  fqipear  to  Inive  a  regular  stiparation.  Through  experimentation,  we  have  found  that 
multiplication  of  the  dettsity  by  yields  an  interesting  result  when  plotted  agaitist  the  transformed  propa¬ 
gation  distatice,  see  Figure  5.  The  vtiriiition  in  atmosirheric  stretigih  between  the  density  curves  is  accoutited 
for  by  this  function  of  ro.  'I’he  combination  of  the  transformation  of  the  propagation  distance  by  jq  and  the 
normalization  by  the  function  appears  to  have  collapsed  the  data  into  predominately  a  single  curve.  The 
results  depicted  in  Figure  5  demonstrate  that  the  branch  point  distribution  is  dependent  on  the  strength  of 
the  turbulence  and  the  distance  of  propagation  wiili  no  other  factors  playing  a  strong  role.  This  implies  that 
the  density  is  a  predictable  parameter  of  the  atmospheric  turbulence  and  should  be  reducahle  to  a  functional 
relationship  similar  to  the  other  paramiUers  used  iti  (luantifying  titmospheric  conditions. 

.\i.  higher  turbulence  strengths  and  greater  proptigatiott  distances  the  density  deviates  from  the  domin.atc 
tretid  suggesting  a  saturation  effect  at  work  in  the  branch  point  disiribtition.  That  the  saturation  effects  appear 
well  below  the  tnaximum  density  that  the  W'FS  could  theoretically  measure  (roughly  GOOOp/. s/m*)  indicates  that 
this  is  some  self  imposed  conslnunt  withiti  the  turbulence. 
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4.5  Empirical  Algebraic  Form  for  Branch  Point  Density 

With  such  a  well  defined  curve  it  is  passible  to  use  tiiis  iiironiiation  to  build  an  empirical  formula  for  the  branch 

1 1/3 

point  density  as  a  function  of  ro  and  propagation  distance,  z.  'I'he  vjdue  of  appears  to  only  be  a 

function  of  the  propagation  distance  beyond  the  branch  point  threshold,  jq-  Examining  the  curve 
versus  s  —  co  in  log-log  space  reveals  a  functional  form  proportional  to 

Incorporating  this  dependence  into  an  empirical  relationshi])  yields 


Stanch  Point  Density  times  rj’^ 


Figure  5.  Branch  point  density  from  single  phase  wheel  data  for  selected  turbulence  strengths.  For  each  curve,  the  turbulence 
stTT.ngth,  given  by  ro,  is  held  constant  while  the  propagation  distance  is  varied  using  the  optical  trombone.  Branch  point 
density  vs  propagation  distance  with  p^,  normalized  by  plotted  versus  scaled  projmgation  distance  (z  —  zo) 


(a)  _  ^  (b) 

Figure  G.  Estimating  the  functional  dependence  of  Pu^r^'^^  on  the  scaled  propagation  distance,  z  —  zq.  (a)  The  density 
data  in  log  —  log  space  shown  as  green  and  blue  dots.  The  green  dots  indicate  those  points  that  appear  to  show  signs  of 
"saturating"  away  from  the  general  behavior  of  the  mam  curves.  Overlaid  with  a  red  line  representing  a  best  linear  fit  that 
demonstrates  the  11/G'**  {lower  dependence,  (b)  The  same  functional  form  overlaid  on  the  data  plot. 


(14) 


PaPr,„,,  —  '^O  ^0  (-  ~  ^0) 


/G 


with  ko  the  wave  number  and  Cp  as  a  scaling  constant.  Also,  an  additional  co  term  was  added  to  account  for 
scaling  and  balancing  the  units  between  the  left  and  right  hand  sides  of  the  equation. 

This  function  is  over-plotted  on  the  measured  branch  point  densities  for  the  data  sets  from  Table  1  in  Figure  7. 
Cp  in  Equation  14  is  set  to  0.7 17,  which  was  determined  throttgh  a  optimization  algorithm  to  fit  each  curve  tising 
a  least  squares  style  norm  and  then  optimizing  the  result  across  all  of  the  curves  simultaneously.  Interestingly, 
this  experimentally  derived  constant  for  Cp  is  almost  exactly  the  inverse  of  the  constant  in  the  Rytov  equation, 
0.5631(2.88/2.91),  or  0.7511.  In  Figure  7  each  subplot  displays  the  density  versus  propagation  distance.  The 
turbulence  strength  increases  from  left  to  right  and  from  top  to  bottom.  The  branch  point  densities  for  each 
configuration  are  plotted  against  the  unsealed  propagation  distance  in  blue.  The  green  curve  represents  what  is 
predicted  from  Equation  14.  The  red  circles  indicate  points  where  the  Rytov  value  was  above  0.4.  As  can  be  seen, 
for  <  0.4,  Equation  14  gives  reasonable  results.  However,  above  «  0.4,  it  fails  and  this  is  approximately 
where  Sasiela'®  notes  that  the  Rytov  approximation  fails  mathematically. 


5.  DISCUSSION 

The  inclusion  of  the  optical  trombone  to  the  standard  ATS  configuration  provided  a  useful  tool  for  the  examina¬ 
tion  of  the  behavior  of  branch  points  within  the  turbulence  generated  by  the  simulator.  The  experiment  provided 
significant  data  into  the  behavior  of  branch  points  within  the  ASALT  lab  optical  systems  and  has  helped  spawn 
a  new  thread  of  research  within  the  labs  regarding  the  aggregate  behavior  of  branch  points  in  our  experiments. 
This  work  is  being  continued  along  several  avenues  and  will  be  the  material  for  a  number  of  future  papers  already 
in  progress. 

Unfortunately,  for  the  original  purpose  of  this  research,  the  discovery  of  saturation  in  the  branch  point 
density  limits  the  use  of  this  parameter  as  a  means  of  verifying  the  simulated  turbulence  from  open  loop  data. 
Further,  that  the  saturation  begins  to  separate  the  measured  densities  from  the  empirically  derived  equation  at 
the  same  Rytov  value,  «  0.4,  that  the  scintillation  index  fails,  the  branch  point  density  alone  offers  no  additional 
calibration  of  the  ATS  than  scintillation  measurements  already  provide. 

Finally  the  work  here  is  valid  only  for  a  single  turbulence  layer  modeled  in  the  .4SALT  laboratory.  Further 
work  is  necessary  to  determine  the  validity  of  these  predictions  in  a  larger  scope. 
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